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We present a theory for disordered interacting electrons that can describe both the Mott and the 
Anderson transition in the respective limits of zero disorder and zero interaction. We use it to 
investigate the T = Mott-Anderson transition at a fixed electron density, as a the disorder strength 
is increased. Surprisingly, we find two critical values of disorder W n fi and W c . For W > W n fl, the 
system enters a "Griffiths" phase, displaying metallic non-Fermi liquid behavior. At even stronger 
disorder, W = W c > W n fi the system undergoes a metal insulator transition, characterized by the 
linear vanishing of both the typical density of states and the typical quasiparticle weight. 



PACS Numbers: 75.20.Hr, 71.55.Jv 

The nature of the metal-insulator transition, is a fun- 
damental problem in condensed matter science. There 
are two basic mechanisms that cause electron localiza- 
tion. Mott demonstrated that electron-electron interac- 
tions, can produce a metal insulator transition (MIT) 
even in a clean system |Q . Anderson discovered that dis- 
order, i. e. strong spatial fluctuations in the potential 
due to impurities Q], can drive a metal insulator transi- 
tion in a system of non interacting electrons. 

Following these early ideas, important advances were 
made following the application of scaling approaches 
[p|-|lO[ to the problem. In the interacting case, these for- 
mulations turned out to be closely connected to Fermi 
liquid ideas [Q. 

These efforts notwithstanding, many basic questions 
remain. In particular, it proved very difficult to incor- 
porate the effects of strong electronic correlations, such 
as the formation of local magnetic moments, in a com- 
prehensive theory of the MIT. This is a serious short- 
coming, since it is well established experimentally that 
the metallic state close to the MIT is characterized by a 
divergent magnetic susceptibility and linear specific heat 
coefficient. These observations form the basis of the two 
fluid phenomenology jnj. 

Very recently, a new approach p2[ to the strong corre- 
lation problem has been developed and successfully ap- 
plied to systems in the vicinity of the Mott transition. 
This dynamical mean-field theory is in its spirit quite 
similar to the well known Bragg- Williams theory of mag- 
netism, and as such becomes exact in the limit of large 
coordination. The approach has furthermore been ex- 
tended to disordered systems [^3|, and used to investi- 
gate phenomena such as disorder-induced local moment 
formation However, if formulated in its strict large- 
coordination limit, the theory misses strong spatial fluc- 
tuations, and thus cannot incorporate Anderson localiza- 



tion effects. 

The goal of the present study is to present a theory that 
can describe both the Mott and the Anderson route to 
localization, and therefore address the interplay of these 
effects. We follow an approach very similar to the well 
known Thouless-Anderson-Palmer (TAP) formulation of 
the mean field theory of spin glasses (l5) . Specifically, we 
treat the correlation aspects of the problem in a dynami- 
cal mean-field theory fashion, but allow spatial variations 
of the order parameter in order to allow for Anderson 
localization effects. The theory is then exact in the non- 
interacting limit, and reduces to the standard dynamical 
mean field theory in absence of disorder. 

For simplicity, we consider a simple single-band Hub- 
bard model with random site energies, as given by the 
Hamiltonian 

H = ^2^2(-tij + £iSij)cl a c jia + U^c\ A c lA c\ lC%A . 

ij a i 

Within the dynamical mean-field theory, all local cor- 
relation functions can be evaluated using a single-site ef- 
fective action of the form 

Sef f (i) = Y,f dT [ dr'c\jT)(S(T - r>) (d T + e* - m) 

Jo Jo 
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Here, we have used functional integration over Grass- 
mann fields Ci i<T (r) that represent electrons of spin a on 
site i, and rii^ir) — cj a {j)ci.a{j)- The " hybridization 
function" Aj(r, t') is obtained by formally integrating 
out all the degrees of freedom on other sites in the lat- 
tice, and is given by 
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The sum over j runs over the z neighbors of the site i, and 
Gj(w n ) =< Cj(ui n )cj(uj n ) > are the local Green's func- 
tions evaluated on site j, but with the site i removed. 
For z finite, and arbitrary lattices, Gj (u) n ) cannot be 
expressed through local Green's functions only, but the 
situation is simpler on a Bethe lattice |l6| ] , where a simple 
recursion relation can be written for this object, express- 
ing it through similar objects on neighboring sites. In 
particular, Ga (u) n ) can be computed from a local action 
of the form identical as in Eq. (2), except that in the 
expression for Aj(r, t'), the sum now runs over z — 1 
neighbors, excluding the site i. 

We note that this local action is identical as the action 
of an Anderson impurity model embedded in a sea of 
conduction electrons described by a hybridization func- 
tion Aj(t, t'). We conclude that the objects G^\uj n ) are 
related by a stochastic recursion relation, that involves 
solving Anderson impurity models with random on-site 
energies £j. 

To make further progress, it is crucial to identify ap- 
propriate order parameters that can characterize differ- 
ent phases of the system and describe quantitatively the 
approach to the transition. In early work, it has already 
been stressed by Anderson Q that a proper description 
of disordered systems should focus on distribution func- 
tions, and that typical rather than the average values 
should be associated with physical observables. Our for- 
malisms maps the original model onto an ensemble of 
Anderson impurity models, and its low energy behavior 
is naturally described in terms of the distribution func- 
tion of the corresponding local density of states (DOS), 
defined as pj — ~ImGj(0) [0. From this distribution 
we can extract the typical DOS pt yp = exp{< hip >}, 
which is a natural order for the metal insulator transi- 
tion. 

On the metallic side of the transition, the distribu- 
tion function of a second quantity, the local quasiparticle 
(QP) weight, which is obtained from the Greens functions 
as qj — -J^jRelGj 1 — Aj]\ u —o, is necessary to character- 
ized the low energy behavior near the transition. Im- 
portant information is obtained from the typical value of 
the random variable <7j, defined as qt yp — exp{< hiqj >}, 
which emerges as a natural order parameter from previ- 
ous studies of the Mott transition. 

It is also useful to consider the average quasiparticle 
(QP) density of states pqp =< Pj/qj > ■ This object is 
very important for thermodynamics, since it is directly 
related to quantities such as the specific heat coefficient 
7 = C/T, or the local spin susceptibility xioc- 

It is instructive to discuss the behavior of these order 
parameters in the previously studied limiting cases. In 
the limit of large lattice coordination spatial fluctuations 
of the bath function Ai(u> n ) are unimportant, and there 
is no qualitative difference between typical and average 
quantities. In the Mott insulating phase there is a gap 



in the density of states, while there is a finite density of 
states on the metallic side of the transition. As the MIT 
is approached from the metallic side, pt yp remains finite, 
but qtyp is found [ p"3| to linearly go to zero. 

Another well studied limit is that of nonintcracting 
electrons on the Bethe lattice, which is known ]2],pj|,p|,^8[ 
to display an Anderson transition. In the Anderson insu- 
lator phase the local density of states has strong spatial 
fluctuations, few sites with discrete bound states near the 
Fermi level have large density of states while the den- 
sity of states in most of the sites is zero. The average 
DOS is finite both in the insulating and in the metallic 
phase, and is non critical at the transition. Similarly, 
by definition q typ = 1 in this noninteracting limit, so it 
also remains non critical. On the other hand, the typical 
density of states pt yp is finite in the metal and zero in 
the Anderson insulator. This quantity is critical, and is 
found to vanishes exponentially |l9) with the distance to 
the transition. 

Equation |^ is a system stochastic equations, i.e. they 
depend on the realization of the random variables de- 
scribing the disorder. To calculate the probability distri- 
butions of pj and qj we use a simulation approach, where 
the probability distribution for the stochastic quantity 
G^ l \u> n ) is sampled from an ensemble of N sites, as orig- 
inally suggested by Abou-Chacra et al. |L6) . To solve An- 
derson impurity models for given bath functions A 7 - (r, r') 
we use the slave boson (SB) mean-field theory Q 
which is known to be qualitatively and even quantita- 
tively correct at low temperature and at low energies. 

We now discuss our results for the nontrivial situation 
where both the disorder and the interactions are present. 
We consider a z = 3 Bethe lattice, in the limit of infinite 
on-site repulsion U at T = and fixed filling n = 0.3, 
in the presence of a uniform distribution of random site 
energies Ei of width W (following the notation of Rcf. 
p6[ , W is measured units of the hopping element t). We 
begin by concentrating on the evolution of the probabil- 
ity distribution of the local quasiparticle weights qi, as 
the disorder is increased. The sites with ^ < 1 rep- 
resent |l4|,[l3| disorder-induced local magnetic moments, 
and as such will dominate the thermodynamic response 
(see the definition of pqp). For weak disorder we expect 
relatively few local moments and the quasiparticle weight 
distribution is peaked at a finite value. As the disorder 
is increased, the distribution of qj-s broadens. At a crit- 
ical value of the disorder W n fh a transition to a NFL 
metallic state takes place. To illustrate this behavior we 
display the integrated distribution of the variable q, n(q) 
for different values of disorder in Fig. 1(a). If n(q) ~ q a , 
as q — > 0, and a < 1, then P(q) — » +oo in this limit. 
Since the local Kondo temperatures T^' ~ qi this 
behavior reflects a singular distribution of Kondo tem- 
peratures. As a result, we immediately obtain non-Fermi 
liquid (NFL) behavior [^2| with diverging 7 and xioc 
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FIG. 1. Evolution of probability distributions for interact- 
ing electrons as a function of disorder at T = 0: (a) integrated 
distribution for local quasiparticle weights (local Kondo tem- 
peratures). Results are presented for W = 1,3,5 (dotted 
lines), W = 7 (dashed line), and W = 9,10,11 (full lines). 
The transition to the NFL regime is signaled by the diver- 
gence of the slope of n(q) at q = 0. (b) The evolution of 
the local DOS distribution is presented by plotting P(ln p) 
for W = 3, 5, 7, 9, 10. We find that the maximum, i. e. 
< In p > shifts, as the transition is approached. Note also 
the extremely large width of the distribution, so that p now 
spans many orders of magnitude. 

at T = 0. As we can see, there is a well denned value 
of disorder W n fi ~ 7, beyond which the slope of n(q) at 
q = diverges, and we enter the NFL phase. It is worth 
mentioning that a similar transition to a NFL metal, well 
before the MIT, has been found from the field-theoretical 
approaches in 2+e dimensions |^-PJl0[|. In the NFL phase 
the thermodynamics is dominated by disorder-induced 
local moments. The probability distribution of the sec- 
ond order parameter p, P(lnp) , for different values of 
the disorder strength is shown in Fig. 1(b). Notice that 
not only the width, but also the maximum of the dis- 
tribution shifts with disorder, a behavior reminiscent of 
an ordinary Anderson transition. The typical DOS is 
strongly depressed at strong disorder. This behavior is 
even more clearly seen if we plot the DOS averages at the 
Fermi energy as a function of disorder, as presented in 
Fig. 2(b). The typical DOS decreases in a clearly linear 
fashion, as the transition at W — W c « 11 is approached. 
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FIG. 2. Order parameters as functions of the disorder 
strength W . In the noninteracting limit (a), the typical DOS 
vanishes exponentially with disorder, while the average DOS 
is non critical. When interactions are present (b), the typical 
DOS decreases linearly with disorder, while at the same time 
the average one diverges. The divergence is clearly seen by 
plotting 1/ < p > a b (dotted line), which vanishes linearly as 
the critical disorder is approached. Both quantities are found 
to be critical at W — W c ~ 11. Also shown is 1/ < p >qp 
(thin full line), which vanishes at W = W n fi ~ 7. Finally, 
we show in (c) the critical behavior of the typical QP weight, 
which also vanishes linearly at W — W c , similarly as in a 
Mott transition 



This should be contrasted |2{| to the U = Anderson 
transition, where we find (see Fig. 2(a)) the decrease to 
be exponential in agreement with analytical results [18|. 
We mention that at least in the noninteracting limit [18|, 
the typical DOS behaves in a fashion which is qualita- 
tively identical to that of the diffusion coefficient. Hav- 
ing this in mind, one is tempted to interpret our results 
as indicating linear behavior of the conductivity near the 
transition, as found experimentally in many "compen- 
sated" systems. Even more dramatic is the behavior of 
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the average DOS which is non-critical both near a con- 
ventional U = Anderson transition, and near a clean 
Mott transition. This quantity is found to diverge at 
the same value of disorder where the typical DOS van- 
ishes. The fact that we indeed have the divergence, is 
further confirmed by plotting 1/ < p > av as a function 
of disorder, as shown by a dotted line in Fig. 2(b). This 
quantity vanishes linearly at the same critical disorder 
W = W c w 11, In the same figure we exhibit the di- 
vergence of the QP DOS, at the transition to the NFL 
phase. Finally, we consider the behavior of qt yp , which 
is also found to vanish linearly at W = W c , similarly as 
in the case of the Mott transition, but in contrast to the 
noninteracting scenario. Physically, this indicates that 
a finite fraction of electrons turn into strictly localized 
magnetic moments at the metal-insulator transition. 

To summarize, in this paper we have presented a new 
self-consistent theory of disordered interacting electrons 
that can describe both the Anderson and the Mott route 
to localization. In this approach, the typical local DOS 
and the typical local resonance width play the a role of 
order parameters, but the entire probability distributions 
are needed to fully characterize the behavior of the sys- 
tem. Our equations take a form of stochastic recursion 
relations for these quantities that involves solving an en- 
semble of Anderson impurity models. As a specific ap- 
plication of this approach, we have considered a large U 
limit of the Hubbard model at a fixed electron density, 
and investigated effects induced by gradually turning on 
the disorder. We find that the correlations effects pro- 
duce dramatic modifications of the conventional Ander- 
son scenario. At intermediate disorder, there is a tran- 
sition to a non-Fermi liquid phase, characterized by sin- 
gular thermodynamics, but conventional transport. At 
larger disorder a metal-insulator transition takes place. 
This is a new type of transition, having some of the fea- 
tures of both the Anderson and the Mott scenario. Re- 
markably, the main features our treatment, a non Fermi 
liquid phase before the metal insulator transition and a 
linearly vanishing conductivity are found in compensated 
doped semiconductors. 

Our framework suggest several research directions. 
One would like to relate response functions that deter- 
mine the transport coefficients to the local order param- 
eters, as was done in the non interacting case by Efetov 
and Viehweger Hj . Our calculations should be extended 
to to the vicinity of half filling where correlations effects 
should be even more pronounced. This study could cast 
some light on the different types of metal insulator tran- 
sitions that occur in compensated and uncompensated 
doped semiconductors. 
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